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Abstract 

In  this  paper  we  present  a class  of  C 2 spatial  interpolating  curves  depending  on  a set 
of  tension  parameters  and  we  illustrate  their  ability  to  reproduce  the  shape  of  the  data. 
The  curves  are  constructed  using  cubic  splines  and  basically  reduce  to  classical  i/-splines 
for  particular  values  of  the  tension  parameters. 


1 Introduction 

Shape-preserving  interpolation  via  functional  as  well  as  parametric  splines  is  a well 
studied  topic  for  the  planar  case.  On  the  other  hand,  shape-preserving  interpolation  for 
spaces  curves  is  considerably  more  complex  than  for  planar  ones  and  the  related  literature 
is  apparently  limited.  On  this  concern,  a considerable  part  of  the  available  schemes  only 
ensures  geometric  continuity  of  the  obtained  curve  (see  [1,  8]  and  references  quoted 
therein).  Recently,  C 2 and  C3  shape- preserving  interpolating  space  curves  have  been 
obtained  using  polynomial  splines  of  variable  degree,  [2,  3,  6].  However,  working  with 
low(fixed)-degree  polynomial  splines  seems  to  be  a standard  choice  in  the  CAD/CAM 
community.  This  motivates  the  careful  investigation  of  shape  preserving  properties  of 
cubic  v— splines  recently  carried  out  in  [7]  and  the  present  paper. 

In  this  paper  we  present  a method  for  constructing  C 2 spatial  interpolating  curves 
reproducing  the  shape  of  the  polygonal  line  which  interpolates  the  given  data.  The  curve 
is  constructed  via  the  so  called  “parametric  approach” , [10] , using  classical  cubic  splines. 
The  shape  of  the  curve  is  controlled  by  the  amplitude  of  the  tangent  vectors  at  the  data 
sites  which  play  the  role  of  tension  parameters.  It  turns  out  that,  for  particular  values 
of  the  tension  parameters,  the  proposed  scheme  provides  a new,  geometrically  evident, 
description  of  classical  C1  — G 2 cubic  splines,  [11],  Moreover,  the  method  produces  a 

suitable  reparameterization  for  the  above  mentioned  curves  ensuring  C2  continuity.  The 
reparameterization  is  a cubic  polynomial  involving  the  tension  parameters  (see  (3.3)). 
Thus,  the  evaluation  of  the  curve  for  a fixed  value  of  the  new  parameter  requires  the 
solution  of  a cubic  equation. 

The  geometric  meaning  of  the  tension  parameters  coupled  with  the  powerful  “shape- 
preserving” properties  of  the  Bernstein-Bezier  representation  can  be  efficiently  used  to 
construct  an  iterative  algorithm  for  C 2 shape-preserving  interpolation.  The  algorithm 
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converges  in  a finite  number  of  iterations  and  requires  at  each  iteration  the  solution  of 
a diagonally  dominant  linear  system. 

The  paper  is  organized  as  follows.  In  Section  2 we  state  the  problem.  In  Section  3 we 
describe  the  construction  of  the  required  interpolant  and  we  illustrate  its  dependence  on 
the  tension  parameters.  The  asymptotic  behavior  and  the  shape-preserving  properties 
of  the  obtained  curve  are  briefly  discussed  in  Section  4.  We  conclude  in  Section  5 with 
a graphical  example. 


2 The  problem 

In  this  section  we  introduce  the  problem  of  shape-preserving  interpolation  by  curves  in 
IR:i . The  adopted  notion  of  shape- preserving  follows  the  definitions  of  [2]  and  [6].  Let 


I i £ 1R3,  i = 0,...,N, 


be  the  interpolation  points  with  I*  / Ii+i . Define,  for  all  admissible  indices, 


L,  :=  L 


+1 


N, 


r Lj_i x 
:=  \ ||L;_1  X 

l 0, 

:=  {ft: 

lo, 


% i 

xLj 

Li 


if  IILi-x  x L 
elsewhere, 

Li+i ' 


>0, 


IILiXL. 


mil 


if  ||Lj_!  X L, 
elsewhere, 


|L i x Lj+i||  > 0, 


where  |a  b c|  denotes  the  determinant  of  the  matrix  with  columns  a,  b,  c.  The  vectors 
N i and  the  scalars  A,  are,  respectively,  the  discrete  binormals  and  the  discrete  torsions 
of  the  data. 

Let  the  parameter  values  cq,  i = 0, . . . , N,  with  cq  < cq+i  be  given,  and  let 


hi  :=  ai+i  - Oi,  i = 0, 1, . . . , N - 1 

be  the  corresponding  spacings.  We  wish  to  construct  a curve  Q(s),  s £ [op, ffw],  which 
interpolates  the  data,  Q (cq)  = I,,  i = 0 such  that  Q £ C2 [<To, (Tn].  In  addi- 
tion, we  also  require  that  Q(s)  is  shape-preserving,  that  is  it  reproduces  the  convexity 
and  torsion  of  the  polygonal  line  connecting  the  interpolation  points.  More  specifically, 
denoting  with  dashes  derivatives  with  respect  to  the  parameter  s,  we  define 


K(s) 


Q '(a)  x Q "(s) 

IIQ'O)!!3 


Q'(s)  ^ 0,  r(s)  := 


Q'(s)  Q"(s)  Qw(s)| 

IIQ'W  x Q"(s)ll2 


if  K(s)  /’  0 


(2.1) 

as  the  curvature  vector  and  the  torsion  of  the  curve  respectively.  Q(s)  is  shape- preserving 
if  it  satisfies  the  following  criteria  ([2,  6,  7]). 


(i)  Convexity  criteria: 

(1.1)  if  Nj  • Ni+i  > 0,  then  K(s)  • N,-  > 0,  j = i,i  + 1,  s £ [cq,cq+i], 

(1.2)  if  Nj  • Ni+i  < 0,  then  K(s)  • Nj,  j — i,i  + 1,  has  one  change  in  sign  in  [cq,<Tj+i], 

(1.3)  if  Ni  • Nj  + 0 then  (K(uj)  • Nj-)(N * • N ,•)  > 0,  j = i - 1,  i,  i + 1. 


(ii)  Torsion  criteria:  if  A,  ^ 0 then  r(s)Aj  > 0,  s £ [cxd- , cri+1] . 


Carla  Manni 


For  the  sake  of  brevity  we  refer  to  [7]  for  the  more  technical  collinearity  and  coplanar- 
ity criteria. 

3 Constructing  the  interpolating  curve 

In  order  to  construct  the  curve  Q we  consider,  as  a first  step,  a cubic  curve  C interpol- 
ating the  data.  We  put 

C(t)ll„i,ai+l]:=Ci(t;X?>\\?)),  (3.1) 

C<(t;  A|o: \ A|1})  :=  I iH$\u)  + Ii+1  H[0)(u)  + (u)  + X^hiT^H^iu), 

t€[(Ti,ai+i],  u:=  (t-  al)/hi,  (3.2) 

where  0 < A-°\a^  < 1 are  shape  parameters,  T,-,  T,+i  are  vectors  to  be  determined 
and  H^\u)  denote  the  elements  of  the  cardinal  basis  for  cubic  Hermite  interpolation, 
that  is  (u)  are  the  polynomials  of  third  degree  such  that 

dlH\3){r) 

dyl  — vljOri'i  ^ 1 1 1* 

One  can  immediately  verify  that  the  curve  (3.2)  interpolates  the  points  I,- , I?+1  at  the 
extremes  of  the  interval  [er^ov+i]  and  has  tangent  vectors  A^T,-,  A^Tj+i  at  the  same 
extremes.  The  parameters  A-°\  A^  determine  the  amplitude  of  the  tangent  vectors  of 
the  curve  at  the  two  end  points  of  the  interval  and  they  control  the  shape  of  the  curve.  To 
be  more  specific,  since  Hq°\u)  + h[°\u)  = 1,  we  have  that  Cj(t;  0, 0)  reduces  to  the  line 
through  Ij,  Ij+i.  Thus,  the  parameters  Af^A^  act  as  tension  parameters  stretching 
the  curve  from  the  classical  Hermite  cubic  interpolating  I,;,  I,+i  with  tangents  T,,  T,;+] 
(A|°\a|^  = 1)  to  the  line  segment  (A-°^,A^  = 0).  The  curve  (3.1)  turns  out  to  be  of 
class  Gl. 

Let  us  consider  now  the  new  global  parameter 

r^Mf.AP)  :=  (TiH^\u)  + ai+lH^\u)  + (3.3) 

X^hiH^iu)  + a{i; )hiHi1l\u). 

It  is  not  difficult  to  see  that,  if 

0<aP),a[1)<1  (3.4) 


dsl(t;Xf'\x<i1)) 


> 0,  t e {oi,ai+i} 


Thus  (3.3)  implicitly  defines  a function  t = t(s),  which  provides  a reparameterization  for 
(3.1).  In  the  following  we  assume  that  conditions  (3.4)  hold  and  we  define 

Q(s)  :=  C(t(s)).  (3.5) 

Since  Q'(erj)  = T,;,  i = 0, . . . , Ar,  Q is  of  class  C1.  For  each  sequence  of  the  tension 
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parameters  A-°\  A,;1'1  we  will  determine  the  tangent  vectors  Tj,Ti+i  so  that  Q is  also 
of  class  C 2 . Let  us  denote  by  dots  derivatives  with  respect  to  the  local  parameter  u. 
Imposing  continuity  of  Q"(s)  at  cq,  i = 1, . . . ,N  — 1,  from  (3.3),  (3.5)  and  from  the 
chain  rule  for  derivatives,  we  obtain 

04%^  - *-i( n/ii-iA^T,  Ci(0+)hjA-o)  — Sj(0+)fe;A-o)T; 

(hi- lA^)3  (hi  Ai0))3 

Thus,  after  some  manipulations,  from  (3.2)  we  have 

+ Tj  + WjTj+i  — Zi,  i = 1, . . . , N — 1,  (3-7) 


h^A^CMf)2 
— ) 

Wi 

hi#' fa- iA^)2 

^2  — )' 

Wi 

Wi  = hj_i(3  - A-°)1)(/iiA-0))2  + /q(3  - x51))(fc4_1A<i>1)a>  (3.8) 

Zi  = —Liihi-^)2  + — Li-ifoAj0*)2. 
uq  Wi 

In  order  to  uniquely  determine  the  vectors  T*  we  need  two  additional  equations  that 
will  be  obtained  by  imposing  boundary  conditions.  Classical  boundary  conditions  are 
periodic  conditions: 

woTjv-i  + To  + woTi  = zo,  «jvTjv_i  + Tjv  + vatTi  = z^v 

(with  Uo,Vo,Un,vn,zo,Zn  defined  according  to  (3.8)  setting  h_i  = h^-i,  A^°]  = A^_i, 

A(i]  = A^_j,  L_i  = Ljv_i,  hN  = h0,  A$  = A(00),  A^  = A^,  Ljy  = L0)  and  end 
tangent  conditions : 

To  = D0,  Tjv  = D;v, 

(where  Do,  D^r  are  given  in  input).  In  the  following  we  will  denote  by  X the  set  of  indices 
{1, . . . , N — 1}  ({0, . . . , N})  when  end  tangent  (periodic)  conditions  are  considered.  It  is 
not  difficult  to  see  that  (3.7)  for  any  choice  of  the  above  mentioned  boundary  conditions 
provide  a diagonally  dominant  system  . 

AT  = z.  (3.9) 

Thus  we  can  state  the  following 

Theorem  3.1  For  any  sequence  A|0\  A,-1^,  i = 0, . . . , IV— 1,  satisfying  (3.4),  there  exists 
a unique  Q £ C2[ao,  &n]  defined  via  (3.1)-(3.3),  (3.5)  which  interpolates  the  given  data 
and  satisfies  periodic  or  end  tangent  conditions. 

We  notice  that  for  A^  = A^  — 1,  system  (3.9)  reduces  to  the  system  for  the 
computation  of  classical  C2  cubic  splines.  Moreover,  if  A^  = A^  = A k,  k £ X,  the 
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curve  C is  of  class  C 1 and  equation  (3.6)  reads 

d2  c / +\  _ 2s’i(0+)  - ) d 

Gi\ai  ) (j^2  )—  \ )• 


dt 2 


dt 


Then  (3.6)  is  equivalent  to  impose  that  the  cubic  curve  (3.1)  is  a Cl-G2  cubic  is-  spline 
[5,  7,  11]  where,  from  (3.3),  for  i £ I 

K2si{ 0+)  - h“21si_j(l-)  (6  - 4Af  - 2A i+i)hf  + (6  - 2A,_,  - 4A Jhr}, 


Vi  := 


(3.10) 


4 Asymptotic  behavior  and  shape-preservation 

In  this  section  we  briefly  discuss  the  asymptotic  behavior  and  the  resulting  shape- 
preserving properties  of  the  curve  Q,  defined  by  (3.1)-(3.3),  (3.5)  and  (3.9),  as  the 
tension  parameters  A.-°\  A-1^  approach  zero.  The  following  lemma  (see  also  [7])  concerns 
the  asymptotic  behavior  of  the  tangents  T,-.  We  omit  the  details  of  the  proof  which  are 
completely  analogous  to  those  of  Theorem  3 in  [9] . 


Lemma  4.1  The  vectors  T *,  i = 0, . . . , N,  obtained  from  (3.9)  are  bounded  independ- 
ently of  Aj°) , Aj1^ , j = 0, . . . , N - 1 . Moreover, 


lim 

A<°>  x<»_>0 


T,  = 


M a|0))2 


L,-i 


(i)  \2 


+ 


hi-i(A._j) 


hi{ Af)2  + hi-1  hi-ri a£\)2  + M A<0)) 


(°)^2  hi 


X bj  — 1 . hi  . 

==:  (1  — Oti)—  h &i~T~ , l £ I. 

Hi- 1 Hi 


(4.1) 


Since  the  tangents  are  bounded  independently  on  the  tension  parameters,  from  the 
previous  section  we  have  that  Q approaches  the  piecewise  linear  function  interpolating 
the  data  as  the  tension  parameters  tend  to  zero.  Moreover,  each  tangent  T*  determined 
by  (3.9)  tends  to  a strictly  convex  combination  of  and  L i/hi  as  the  tension 

parameters  A^°\,  A^  tend  to  zero  while  A-i\/A.0)  remains  bounded  and  strictly  positive. 
Due  to  these  two  main  facts,  we  are  able  to  easy  control  the  shape  of  the  curve  Q and 
to  ensure  that  it  reproduces  the  shape  of  the  data  as  the  tension  parameters  approach 
zero  as  we  will  discuss  briefly  in  the  following. 

Since  C and  Q only  differ  for  a reparameterization  they  have  the  same  image.  Thus,  as 
far  as  the  shape-preserving  properties  are  concerned,  we  can  consider  the  expression  of  C. 
As  noticed  in  Section  3,  if  Aji\  = A[0),  i £ I,  the  curve  C with  Tj  obtained  by  (3.9),  is  a 
Cl-G 2 cubic  ^-spline.  In  such  a case,  using  (3.10),  the  careful  shape  analysis  carried  out  in 
[7]  and  the  resulting  algorithm  can  be  considered.  However,  the  simple  geometric  meaning 
of  the  tension  parameters  A-1'  coupled  with  the  “shape-preserving”  properties  of  the 
Bezier-Bernstein  representation,  allow  us  to  more  easily  establish  the  shape-preserving 
results  also  for  completely  general  configurations  of  A-l\,  A^°\  Thus,  we  express  the 
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curve  segment  Ci(t;  x\°\  X^)  in  Bezier-Bernstein  form: 

3 


C«(t;  Af >,  A«)  = (f )tl(l  - t)3~l, 


1=0 


Cifi  :=  If,  Citi  :=  I * + IhiX^Ti,  C,l2  :=  Ii+i  - §M4(1)Ti+1,  Ci)3  :=  Ii+i. 
Let  us  consider  at  the  beginning  the  convexity  criteria. 


Lemma  4.2  If  Nt  • N ,•  ^ 0 and 


(i) 


A<°> 


c > 0,  then 


lim  (K(o-j)  • Nj)(Ni  • Nj)  > 0. 


Proof:  From  the  properties  of  Bezier  curves  (see  [5])  and  from  (2.1)  and  (3.5) 
sgn(K (cTj)  • Nj)  = sgn((Ci,i  - Cij0)  x (Ci)2  - C^i))  • Nj 


sgn 


_ /_  A \0)hir_  X^h;  ' 

Tj  x ^ Tj  Tj+i 


•No- 


where sgn (y)  denotes  the  sign  of  y.  Moreover,  from  (4.1) 


lim  (T i x Lj)  • Nj  = 


a*  77  x Lj  + (1  - a*)^— ^ x L * 

hi  hi- 1 


N,  = (1,  Qi)N,: 


hi  — 


N,. 


Hence,  we  obtain  the  assertion  if  N*  • Nj  / 0.  □ 

The  previous  lemma  ensures  that,  if  A-°\,  A-1^  are  small  enough  the  third  convexity 
criterion,  (i.3),  stated  in  Section  2 is  satisfied.  In  addition,  the  sign  of  K(cXfc)-Nj,  k = i,i+ 
1 can  be  checked  considering  the  Bezier  coefficients  C iti,  l = 0,1, 2, 3,  of  C*.  Furthermore, 
thanks  to  the  shape-preserving  properties  of  totally  positive  bases,  for  small  values  of  the 
tension  parameters,  (see  [4])  the  number  of  changes  in  sign  of  K(s)  • Nj,  s 6 [oj,  oj+i] 
is  bounded  by  the  number  of  changes  of  sign  in  the  pair  K(ofc)  ■ Nj,k  = i,i  + 1.  Thus, 
also  the  first  and  the  second  convexity  criteria  (i.l)  and  (i.2)  are  satisfied  if  the  tension 
parameters  are  small  enough. 

As  far  as  the  torsion  is  concerned,  we  recall  that  the  sign  of  the  torsion  of  a cubic 
curve  coincides  with  the  sign  of  the  discrete  torsion  of  its  Bezier  control  polygon  (see  for 
example  [5])  thus  it  is  not  difficult  to  obtain  the  following 

Ae-> 

Lemma  4.3  If  Aj  ^ 0 and  —>  c > 0,  j = i,  i + 1,  then 

Xj 


r(s)Aj  > 0,  s e [cr+, 


With  similar  arguments  it  is  not  difficult  to  prove  that  also  the  collinearity  and  the 
coplanarity  criteria  stated  in  [7]  are  fulfilled  as  the  tension  parameters  approach  zero. 
We  omit  the  details  for  the  sake  of  brevity. 

Summarizing,  from  the  previous  discussion  it  follows  that  if  the  tension  parameters  are 
small  enough  then  the  Bezier  control  polygon  of  C reproduces  the  shape  of  the  data  and 
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the  curve  C does  the  same  thanks  to  the  properties  of  Bezier-Bernstein  representation. 
Thus,  to  obtain  an  automatic  algorithm  to  compute  the  C2  interpolant  Q defined  by 
(3.5),  satisfying  convexity  and  torsion  criteria,  basically  we  have  to  perform  the  following 
steps: 

(a)  for  a given  sequence  of  the  tension  parameters  solve  the  system  (3.9)  and  compute 
the  Bezier  coefficients  of  the  resulting  curve  C; 

(b)  check  if  the  control  polygon  of  each  segment  C?;  satisfies  the  convexity  and  torsion 
criteria; 

(c)  if  this  is  not  the  case  reduce  the  values  of  the  related  tension  parameters  according 
to  a given  rule  and  go  to  step  (a). 

5 A graphical  example 

To  illustrate  the  performance  of  the  presented  scheme  we  consider  the  data  proposed 
in  [7],  Example  2,  consisting  of  20  points  with  uniform  parameterization  in  [0,1].  End 
tangent  boundary  conditions  have  been  used  (see  Table  2 in  [7]).  Figures  1-3  show  the 
behavior  of  the  obtained  C 2 curve  Q compared  with  the  classical  C2  cubic  spline.  The 
shape-preserving  curve  Q is  defined  by  the  following  sequence  of  tension  parameters 

A-0)  : .6  .6  1 .9  .9  1 1 1 1 1 1 1 1 1 .75  1 1 1 1 

a|1}  : .9  .6  .6  1 .9  .9  1 1 1 1 1 1 1 1 1 1 1 .75  1. 


FlG.  1.  C 2 cubic  spline  (left)  and  Q (right). 
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